Single Cell Atlas: a single-cell multi-omics human cell encyclopedia.

Single-cell sequencing datasets are key in biology and medicine for unraveling insights into heterogeneous cell populations with unprecedented resolution. Here, we construct a single-cell multi-omics map of human tissues through in-depth characterizations of datasets from five single-cell omics, spatial transcriptomics, and two bulk omics across 125 healthy adult and fetal tissues. We construct its complement web-based platform, the Single Cell Atlas (SCA, www.singlecellatlas.org ), to enable vast interactive data exploration of deep multi-omics signatures across human fetal and adult tissues. The atlas resources and database queries aspire to serve as a one-stop, comprehensive, and time-effective resource for various omics studies.

scDrugPrio: a framework for the analysis of single-cell transcriptomics to address multiple problems in precision medicine in immune-mediated inflammatory diseases.

BACKGROUND: Ineffective drug treatment is a major problem for many patients with immune-mediated inflammatory diseases (IMIDs). Important reasons are the lack of systematic solutions for drug prioritisation and repurposing based on characterisation of the complex and heterogeneous cellular and molecular changes in IMIDs. METHODS: Here, we propose a computational framework, scDrugPrio, which constructs network models of inflammatory disease based on single-cell RNA sequencing (scRNA-seq) data. scDrugPrio constructs detailed network models of inflammatory diseases that integrate information on cell type-specific expression changes, altered cellular crosstalk and pharmacological properties for the selection and ranking of thousands of drugs. RESULTS: scDrugPrio was developed using a mouse model of antigen-induced arthritis and validated by improved precision/recall for approved drugs, as well as extensive in vitro, in vivo, and in silico studies of drugs that were predicted, but not approved, for the studied diseases. Next, scDrugPrio was applied to multiple sclerosis, Crohn's disease, and psoriatic arthritis, further supporting scDrugPrio through prioritisation of relevant and approved drugs. However, in contrast to the mouse model of arthritis, great interindividual cellular and gene expression differences were found in patients with the same diagnosis. Such differences could explain why some patients did or did not respond to treatment. This explanation was supported by the application of scDrugPrio to scRNA-seq data from eleven individual Crohn's disease patients. The analysis showed great variations in drug predictions between patients, for example, assigning a high rank to anti-TNF treatment in a responder and a low rank in a nonresponder to that treatment. CONCLUSIONS: We propose a computational framework, scDrugPrio, for drug prioritisation based on scRNA-seq of IMID disease. Application to individual patients indicates scDrugPrio's potential for personalised network-based drug screening on cellulome-, genome-, and drugome-wide scales. For this purpose, we made scDrugPrio into an easy-to-use R package ( https://github.com/SDTC-CPMed/scDrugPrio ).

An interactive atlas of genomic, proteomic, and metabolomic biomarkers promotes the potential of proteins to predict complex diseases.

Multiomics analyses have identified multiple potential biomarkers of the incidence and prevalence of complex diseases. However, it is not known which type of biomarker is optimal for clinical purposes. Here, we make a systematic comparison of 90 million genetic variants, 1,453 proteins, and 325 metabolites from 500,000 individuals with complex diseases from the UK Biobank. A machine learning pipeline consisting of data cleaning, data imputation, feature selection, and model training using cross-validation and comparison of the results on holdout test sets showed that proteins were most predictive, followed by metabolites, and genetic variants. Only five proteins per disease resulted in median (min-max) areas under the receiver operating characteristic curves for incidence of 0.79 (0.65-0.86) and 0.84 (0.70-0.91) for prevalence. In summary, our work suggests the potential of predicting complex diseases based on a limited number of proteins. We provide an interactive atlas (macd.shinyapps.io/ShinyApp/) to find genomic, proteomic, or metabolomic biomarkers for different complex diseases.

Cardiovascular Consequences of Uremic Metabolites: an Overview of the Involved Signaling Pathways.

The crosstalk of the heart with distant organs such as the lung, liver, gut, and kidney has been intensively approached lately. The kidney is involved in (1) the production of systemic relevant products, such as renin, as part of the most essential vasoregulatory system of the human body, and (2) in the clearance of metabolites with systemic and organ effects. Metabolic residue accumulation during kidney dysfunction is known to determine cardiovascular pathologies such as endothelial activation/dysfunction, atherosclerosis, cardiomyocyte apoptosis, cardiac fibrosis, and vascular and valvular calcification, leading to hypertension, arrhythmias, myocardial infarction, and cardiomyopathies. However, this review offers an overview of the uremic metabolites and details their signaling pathways involved in cardiorenal syndrome and the development of heart failure. A holistic view of the metabolites, but more importantly, an exhaustive crosstalk of their known signaling pathways, is important for depicting new therapeutic strategies in the cardiovascular field.

Reversal of Pulmonary Hypertension in a Human-Like Model: Therapeutic Targeting of Endothelial DHFR.

BACKGROUND: Pulmonary hypertension (PH) is a progressive disorder characterized by remodeling of the pulmonary vasculature and elevated mean pulmonary arterial pressure, resulting in right heart failure. METHODS: Here, we show that direct targeting of the endothelium to uncouple eNOS (endothelial nitric oxide synthase) with DAHP (2,4-diamino 6-hydroxypyrimidine; an inhibitor of GTP cyclohydrolase 1, the rate-limiting synthetic enzyme for the critical eNOS cofactor tetrahydrobiopterin) induces human-like, time-dependent progression of PH phenotypes in mice. RESULTS: Critical phenotypic features include progressive elevation in mean pulmonary arterial pressure, right ventricular systolic blood pressure, and right ventricle (RV)/left ventricle plus septum (LV+S) weight ratio; extensive vascular remodeling of pulmonary arterioles with increased medial thickness/perivascular collagen deposition and increased expression of PCNA (proliferative cell nuclear antigen) and alpha-actin; markedly increased total and mitochondrial superoxide production, substantially reduced tetrahydrobiopterin and nitric oxide bioavailabilities; and formation of an array of human-like vascular lesions. Intriguingly, novel in-house generated endothelial-specific dihydrofolate reductase (DHFR) transgenic mice (tg-EC-DHFR) were completely protected from the pathophysiological and molecular features of PH upon DAHP treatment or hypoxia exposure. Furthermore, DHFR overexpression with a pCMV-DHFR plasmid transfection in mice after initiation of DAHP treatment completely reversed PH phenotypes. DHFR knockout mice spontaneously developed PH at baseline and had no additional deterioration in response to hypoxia, indicating an intrinsic role of DHFR deficiency in causing PH. RNA-sequencing experiments indicated great similarity in gene regulation profiles between the DAHP model and human patients with PH. CONCLUSIONS: Taken together, these results establish a novel human-like murine model of PH that has long been lacking in the field, which can be broadly used for future mechanistic and translational studies. These data also indicate that targeting endothelial DHFR deficiency represents a novel and robust therapeutic strategy for the treatment of PH.

Comprehensive multiscale analysis of lactate metabolic dynamics in vitro and in vivo using highly responsive biosensors.

As a key glycolytic metabolite, lactate has a central role in diverse physiological and pathological processes. However, comprehensive multiscale analysis of lactate metabolic dynamics in vitro and in vivo has remained an unsolved problem until now owing to the lack of a high-performance tool. We recently developed a series of genetically encoded fluorescent sensors for lactate, named FiLa, which illuminate lactate metabolism in cells, subcellular organelles, animals, and human serum and urine. In this protocol, we first describe the FiLa sensor-based strategies for real-time subcellular bioenergetic flux analysis by profiling the lactate metabolic response to different nutritional and pharmacological conditions, which provides a systematic-level view of cellular metabolic function at the subcellular scale for the first time. We also report detailed procedures for imaging lactate dynamics in live mice through a cell microcapsule system or recombinant adeno-associated virus and for the rapid and simple assay of lactate in human body fluids. This comprehensive multiscale metabolic analysis strategy may also be applied to other metabolite biosensors using various analytic platforms, further expanding its usability. The protocol is suited for users with expertise in biochemistry, molecular biology and cell biology. Typically, the preparation of FiLa-expressing cells or mice takes 2 days to 4 weeks, and live-cell and in vivo imaging can be performed within 1-2 hours. For the FiLa-based assay of body fluids, the whole measuring procedure generally takes ~1 min for one sample in a manual assay or ~3 min for 96 samples in an automatic microplate assay.

Structurally-informed human interactome reveals proteome-wide perturbations by disease mutations.

Human genome sequencing studies have identified numerous loci associated with complex diseases. However, translating human genetic and genomic findings to disease pathobiology and therapeutic discovery remains a major challenge at multiscale interactome network levels. Here, we present a deep-learning-based ensemble framework, termed PIONEER (Protein-protein InteractiOn iNtErfacE pRediction), that accurately predicts protein binding partner-specific interfaces for all known protein interactions in humans and seven other common model organisms, generating comprehensive structurally-informed protein interactomes. We demonstrate that PIONEER outperforms existing state-of-the-art methods. We further systematically validated PIONEER predictions experimentally through generating 2,395 mutations and testing their impact on 6,754 mutation-interaction pairs, confirming the high quality and validity of PIONEER predictions. We show that disease-associated mutations are enriched in PIONEER-predicted protein-protein interfaces after mapping mutations from ~60,000 germline exomes and ~36,000 somatic genomes. We identify 586 significant protein-protein interactions (PPIs) enriched with PIONEER-predicted interface somatic mutations (termed oncoPPIs) from pan-cancer analysis of ~11,000 tumor whole-exomes across 33 cancer types. We show that PIONEER-predicted oncoPPIs are significantly associated with patient survival and drug responses from both cancer cell lines and patient-derived xenograft mouse models. We identify a landscape of PPI-perturbing tumor alleles upon ubiquitination by E3 ligases, and we experimentally validate the tumorigenic KEAP1-NRF2 interface mutation p.Thr80Lys in non-small cell lung cancer. We show that PIONEER-predicted PPI-perturbing alleles alter protein abundance and correlates with drug responses and patient survival in colon and uterine cancers as demonstrated by proteogenomic data from the National Cancer Institute's Clinical Proteomic Tumor Analysis Consortium. PIONEER, implemented as both a web server platform and a software package, identifies functional consequences of disease-associated alleles and offers a deep learning tool for precision medicine at multiscale interactome network levels.

Machine learning-based integration of network features and chemical structure of compounds for SARS-CoV-2 drug effect analysis.

High drug development costs and the limited number of new annual drug approvals increase the need for innovative approaches for drug effect prediction. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the cause of coronavirus disease 2019 (COVID-19), led to a global pandemic with high morbidity and mortality. Although effective preventive measures exist, there are few effective treatments for hospitalized patients with SARS-CoV-2 infection. Drug repurposing and drug effect prediction are promising strategies that could shorten development time and reduce costs compared with de novo drug discovery. In this work, we present a machine learning framework to integrate a variety of target network features and physicochemical properties of compounds, and analyze their influence on the therapeutic effects for SARS-CoV-2 infection and on host cell cytotoxic effects. Random forest models trained on compounds with known experimental effects on SARS-CoV-2 infection and subsequent feature importance analysis based on Shapley values provided insights into the determinants of drug efficacy and cytotoxicity, which can be incorporated into novel drug discovery approaches. Given the complexity of molecular mechanisms of drug action and limited sample sizes, our models achieve a reasonable mean area under the receiver operating characteristic curve (ROC-AUC) of 0.73 on an unseen validation set. To our knowledge, this is the first work to incorporate a combination of network and physicochemical features of compounds into a machine learning model to predict drug effects on SARS-CoV-2 infection. Our systems pharmacology-based machine learning framework can be used to classify other existing drugs for SARS-CoV-2 infection and can easily be adapted to drug effect prediction for future viral outbreaks.

Cord blood DNA methylation signatures associated with preeclampsia are enriched for cardiovascular pathways: insights from the VDAART trial.

BACKGROUND: Preeclampsia has been associated with maternal epigenetic changes, in particular DNA methylation changes in the placenta. It has been suggested that preeclampsia could also cause DNA methylation changes in the neonate. We examined DNA methylation in relation to gene expression in the cord blood of offspring born to mothers with preeclampsia. METHODS: This study included 128 mother-child pairs who participated in the Vitamin D Antenatal Asthma Reduction Trial (VDAART), where assessment of preeclampsia served as secondary outcome. We performed an epigenome-wide association study of preeclampsia and cord blood DNA methylation (Illumina 450 K chip). We then examined gene expression of the same subjects for validation and replicated the gene signatures in independent DNA methylation datasets. Lastly, we applied functional enrichment and network analyses to identify biological pathways that could potentially be involved in preeclampsia. FINDINGS: In the cord blood samples (n = 128), 263 CpGs were differentially methylated (FDR <0.10) in preeclampsia (n = 16), of which 217 were annotated. Top pathways in the functional enrichment analysis included apelin signaling pathway and other endothelial and cardiovascular pathways. Of the 217 genes, 13 showed differential expression (p's < 0.001) in preeclampsia and 11 had been previously related to preeclampsia (p's < 0.0001). These genes were linked to apelin, cGMP and Notch signaling pathways, all having a role in angiogenic process and cardiovascular function. INTERPRETATION: Preeclampsia is related to differential cord blood DNA methylation signatures of cardiovascular pathways, including the apelin signaling pathway. The association of these cord blood DNA methylation signatures with offspring's long-term morbidities due to preeclampsia should be further investigated. FUNDING: VDAART is funded by National Heart, Lung, and Blood Institute grants of R01HL091528 and UH3OD023268. HMK is supported by Jane and Aatos Erkko Foundation, Paulo Foundation, and the Pediatric Research Foundation. HM is supported by K01 award from NHLBI (1K01HL146977-01A1). PK is supported by K99HL159234 from NIH/NHLBI.

scDrugPrio: A framework for the analysis of single-cell transcriptomics to address multiple problems in precision medicine in immune-mediated inflammatory diseases.

Background: Ineffective drug treatment is a major problem for many patients with immune-mediated inflammatory diseases (IMIDs). Important reasons are the lack of systematic solutions for drug prioritisation and repurposing based on characterisation of the complex and heterogeneous cellular and molecular changes in IMIDs. Methods: Here, we propose a computational framework, scDrugPrio, which constructs network models of inflammatory disease based on single-cell RNA sequencing (scRNA-seq) data. scDrugPrio constructs detailed network models of inflammatory diseases that integrate information on cell type-specific expression changes, altered cellular crosstalk and pharmacological properties for the selection and ranking of thousands of drugs. Results: scDrugPrio was developed using a mouse model of antigen-induced arthritis and validated by improved precision/recall for approved drugs, as well as extensive in vitro, in vivo, and in silico studies of drugs that were predicted, but not approved, for the studied diseases. Next, scDrugPrio was applied to multiple sclerosis, Crohn's disease, and psoriatic arthritis, further supporting scDrugPrio through prioritisation of relevant and approved drugs. However, in contrast to the mouse model of arthritis, great interindividual cellular and gene expression differences were found in patients with the same diagnosis. Such differences could explain why some patients did or did not respond to treatment. This explanation was supported by the application of scDrugPrio to scRNA-seq data from eleven individual Crohn's disease patients. The analysis showed great variations in drug predictions between patients, for example, assigning a high rank to anti-TNF treatment in a responder and a low rank in a nonresponder to that treatment. Conclusion: We propose a computational framework, scDrugPrio, for drug prioritisation based on scRNA-seq of IMID disease. Application to individual patients indicates scDrugPrio's potential for personalised network-based drug screening on cellulome-, genome-, and drugome-wide scales. For this purpose, we made scDrugPrio into an easy-to-use R package (https://github.com/SDTC-CPMed/scDrugPrio).

Integration of circulating microRNAs and transcriptome signatures identifies early-pregnancy biomarkers of preeclampsia.

BACKGROUND: MicroRNAs (miRNAs) have been implicated in the pathobiology of preeclampsia, a common hypertensive disorder of pregnancy. In a nested matched case-control cohort within the Vitamin D Antenatal Asthma Reduction Trial (VDAART), we previously identified peripheral blood mRNA signatures related to preeclampsia and vitamin D status (≤30 ng/mL) during gestation from 10 to 18 weeks, using differential expression analysis. METHODS: Using quantitative PCR arrays, we conducted profiling of circulating miRNAs at 10-18 weeks of gestation in the same VDAART cohort to identify differentially expressed (DE) miRNAs associated with preeclampsia and vitamin D status. For the validation of the expression of circulating miRNA signatures in the placenta, the HTR-8/SVneo trophoblast cell line was used. Targets of circulating miRNA signatures in the preeclampsia mRNA signatures were identified by consensus ranking of miRNA-target prediction scores from four sources. The connected component of target signatures was identified by mapping to the protein-protein interaction (PPI) network and hub targets were determined. As experimental validation, we examined the gene and protein expression of IGF1R, one of the key hub genes, as a target of the DE miRNA, miR-182-5p, in response to a miR-182-5p mimic in HTR-8/SVneo cells. RESULTS: Pregnant women with preeclampsia had 16 circulating DE miRNAs relative to normal pregnancy controls that were also DE under vitamin D insufficiency (9/16 = 56% upregulated, FDR < .05). Thirteen miRNAs (13/16 = 81.3%) were detected in HTR-8/SVneo cells. Overall, 16 DE miRNAs had 122 targets, of which 87 were unique. Network analysis demonstrated that the 32 targets of DE miRNA signatures created a connected subnetwork in the preeclampsia module with CXCL8, CXCL10, CD274, MMP9 and IGF1R having the highest connectivity and centrality degree. In an in vitro validation experiment, the introduction of an hsa-miR-182-5p mimic resulted in significant reduction of its target IGF1R gene and protein expression within HTR-8/SVneo cells. CONCLUSIONS: The integration of the circulating DE miRNA and mRNA signatures associated preeclampsia added additional insights into the subclinical molecular signature of preeclampsia. Our systems and network biology approach revealed several biological pathways, including IGF-1, that may play a role in the early pathophysiology of preeclampsia. These pathways and signatures also denote potential biomarkers for the early stages of preeclampsia and suggest possible preventive measures.

MYC overrides HIF-1α to regulate proliferating primary cell metabolism in hypoxia.

Hypoxia requires metabolic adaptations to sustain energetically demanding cellular activities. While the metabolic consequences of hypoxia have been studied extensively in cancer cell models, comparatively little is known about how primary cell metabolism responds to hypoxia. Thus, we developed metabolic flux models for human lung fibroblast and pulmonary artery smooth muscle cells proliferating in hypoxia. Unexpectedly, we found that hypoxia decreased glycolysis despite activation of hypoxia-inducible factor 1α (HIF-1α) and increased glycolytic enzyme expression. While HIF-1α activation in normoxia by prolyl hydroxylase (PHD) inhibition did increase glycolysis, hypoxia blocked this effect. Multi-omic profiling revealed distinct molecular responses to hypoxia and PHD inhibition, and suggested a critical role for MYC in modulating HIF-1α responses to hypoxia. Consistent with this hypothesis, MYC knockdown in hypoxia increased glycolysis and MYC over-expression in normoxia decreased glycolysis stimulated by PHD inhibition. These data suggest that MYC signaling in hypoxia uncouples an increase in HIF-dependent glycolytic gene transcription from glycolytic flux.

Drugst.One - A plug-and-play solution for online systems medicine and network-based drug repurposing.

In recent decades, the development of new drugs has become increasingly expensive and inefficient, and the molecular mechanisms of most pharmaceuticals remain poorly understood. In response, computational systems and network medicine tools have emerged to identify potential drug repurposing candidates. However, these tools often require complex installation and lack intuitive visual network mining capabilities. To tackle these challenges, we introduce Drugst.One, a platform that assists specialized computational medicine tools in becoming user-friendly, web-based utilities for drug repurposing. With just three lines of code, Drugst.One turns any systems biology software into an interactive web tool for modeling and analyzing complex protein-drug-disease networks. Demonstrating its broad adaptability, Drugst.One has been successfully integrated with 21 computational systems medicine tools. Available at https://drugst.one, Drugst.One has significant potential for streamlining the drug discovery process, allowing researchers to focus on essential aspects of pharmaceutical treatment research.